Noise correlations in a flux qubit with tunable tunnel coupling 
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We have measured flux- noise correlations in a tunable superconducting flux qubit. The device 
consists of two loops that independently control the qubit's energy splitting and tunnel coupling. 
Low frequency flux noise in the loops causes fluctuations of the qubit frequency and leads to de- 
phasing. Since the noises in the two loops couple to different terms of the qubit Hamiltonian, a 
measurement of the dephasing rate at different bias points provides a way to extract both the am- 
plitude and the sign of the noise correlations. We find that the flux fluctuations in the two loops 
are anti-correlated, consistent with a model where the flux noise is generated by randomly oriented 
unpaired spins on the metal surface. 



The performance of qubits based on superconducting 
circuits is affected by the ubiquitous low-frequency flux 
noise present in such devices. Early measurements on 
superconducting quantum interference devices (SQUIDs) 
showed the existence of flux noise with a l//-like power 
spectrum [1 . The noise was found to be universal in the 
sense that it only depended weakly on materials and sam- 
ple dimensions [2 . In more recent years, it has become 
clear that the same excess flux noise limits the phase co- 
herence of superconducting qubits. A slowly varying flux 
leads to fluctuations of the qubit energy levels, which has 
been shown to be the dominant dephasing mechanism in 
both flux [3H5] and phase qubits [6 . Finding the origin 
and a possible remedy for the excess flux noise is there- 
fore of importance to facilitate further improvements of 
intrinsic coherence times in these devices. 

Theories have been developed to explain the micro- 
scopic origins of the noise, involving the existence of real 
or effective unpaired spins in the vicinity of the super- 
conducting structures [THIQ]. Various experiments have 
shown the presence of a large number of localized spins 
on the surface of thin films of normal [TT and supercon- 
ducting metals [12j, as well as in Si/Si02 interfaces [13 . 
Sendelbach et al. measured correlations between flux and 
inductance noise in SQUIDs, possibly suggesting the for- 
mation of spin clusters [14 . The results suggest that the 
flux noise is related to the non-equilibrium dynamics of 
the spin system, possibly described by spin glass models 
[15] or fractal spin clusters [16]. Recent measurements 
on qubits with different geometries indicate that the flux 
noise scales as (l/w)^ in agreement with the surface spin 
model [17] (/ is the length and w is the width of the 
superconducting wires). 

Flux noise correlations have been studied in a system 
of coupled qubits sharing parts of their loops [18]. It 
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FIG. 1: (a) A standard flux qubit. The two diabatic states 
correspond to clockwise and counter-clockwise circulating cur- 
rents, (b) A tunable flux qubit. The smallest junction has 
been replaced by an additional loop, (c) The actual sample 
design. The qubit is embedded in a SQUID used for read-out 
of the qubit state. The line widths are not drawn to scale. 



was found that the flux fluctuations originating from the 
shared branch lead to correlations in the noise of the two 
qubits. In this work, we use a single, two-loop qubit to 
investigate flux noise correlations between different parts 
within a single qubit [Fig. [ijb)]. The fluxes in the two 
loops couple to the longitudinal and transverse compo- 
nents of the qubit's Hamiltonian, respectively. The abil- 
ity to control both parameters of the Hamiltonian allows 
us to investigate qubit coherence properties at different 
frequencies while remaining at the optimal bias point 
[19 . By measuring the qubit dephasing rate as a func- 
tion of flux bias, we can extract the correlations between 
the flux fluctuations in the two loops. Knowing that the 
two fluxes couple differently to the qubit energy allows 
us to determine the sign of the correlations. We find the 
flux noise to be strongly anti-correlated, in agreement 
with a model where the noise is generated by spins on 
the superconductor surface. 

The standard flux qubit consists of a superconducting 
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loop with three or more Josephson junctions [Fig. [TJa)]. 
The diabatic states correspond to clockwise and coun- 
terclockwise circulating currents, respectively [20l [21]. 
When the flux in the loop is close to half a flux quantum, 
$ = 0.5^0 7 the Hamiltonian is H = — (^crx + |crz) within 
a two- level approximation. Here, e = 2/p(^ — 0.5$o) 
where Ip is the persistent current and = h/2e. The 
tunnel coupling A is fixed by fabrication and determined 
by the size of the smallest Josephson junction. By replac- 
ing the smallest junction with a second loop, the tunnel 
coupling A can be tuned in situ by the flux ^2 = f2^o 
in the second loop [Fig. [TJb)] |T9l|22]. Due to the device 
geometry, the detuning e is controlled by the combined 
flux (/i -\- f 2/2)^0^ where fi = ^i/^o is the normalized 
flux in loop 1 [21 . To clarify the relation between the 
qubit parameters and the two fluxes, we introduce the 
effective fluxes 

/z = (/l+/2/2)-l/2, A = /2. (1) 

The energy separation of the qubit is then given by 
Eoi/h = T^^Ta^, with e = A = A(/,). (2) 

The qubit is embedded in a dc SQUID for reading out the 
qubit state. The read-out is implemented by applying a 
short current pulse to the SQUID; due to the inductive 
coupling between the SQUID and the qubit, the SQUID 
switching probability will vary depending on the qubit 
state [23]. 

The design of the actual device is shown in Fig. [TJc) 
(see the supplementary material for a detailed drawing of 
the device [24]). The structure is made of aluminum and 
is fabricated using e-beam lithography and shadow evap- 
oration. There are two local current-bias lines, which, 
together with the global external field, allows us to con- 
trol the fluxes in the two qubit loops and in the read- 
out SQUID independently. Due to the close proximity 
of the loops and the bias lines there is substantial cross- 
coupling between the different elements. During the mea- 
surements we applied appropriate compensation currents 
to the bias lines to compensate the unwanted couplings. 
The residual unwanted coupling is less than 1%. 

Figure [2ja) shows two qubit spectra measured for two 
values of /x- The measurement was done by applying a 
long microwave pulse (3 jas) to saturate the qubit before 
reading out its state. The spectra have the form expected 
from Eq. but with different values of A. The horizon- 
tal feature around 3.2 GHz is due to the plasma mode of 
the read-out SQUID. In Fig. [2|b) we plot the measured 
A as a function of /x. The tunnel coupling A is sym- 
metric around /x = and can be tuned experimentally 
from 1.5 to 15 GHz. The dashed line shows the results 
of a numerical simulation of the qubit energy levels for 
realistic fabrication parameters. To quantify the relation 
between A and /x we plot the numerical derivative of 
the data for a region around /x = 0.39 [Fig. [2];c)]. We 
have marked the positions of the two operating points 
with A = 2.5 GHz and A = 3.5 GHz where we will per- 
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FIG. 2: (a) Qubit spectra, measured for two different val- 
ues of /x. The features at 3.2 GHz are due to the SQUID 
plasma frequency, (b) Qubit tunnel coupling A measured as 
a function of /x. The dashed line is the result of a numerical 
simulation, (c) Numerical derivative dA/df^ of the data in 

(b). 



form coherence measurements. The sensitivity dA/df^ 
is different for those two points. 

The results of Fig. [2] show that our device behaves 
as expected and that we can tune the qubit parameters 
A and e independently by applying fluxes in the two 
loops. To characterize the qubit coherence as a func- 
tion of those parameters, we perform energy-relaxation 
(Ti), free-induction decay (FID) and Hahn-echo mea- 
surements. Ti is measured by applying a tt pulse to the 
qubit and delaying the read-out. Figure |3|a) shows the 
extracted Ti decay times as a function of detuning /z, 
measured for the two values of A marked in Fig.|2|c). Ti 
is fairly independent of /z. The difference in Ti for the 
two values of A can be attributed to differences in detun- 
ing from the SQUID plasma mode at /piasma ~ 3. 2 GHz. 

The FID sequence consists of two 7r/2 pulses separated 
by a time t. The FID is sensitive to low-frequency noise 
that causes adiabatic fluctuations in the qubit energy 
splitting £^01 • The Hahn-echo contains an extra tt pulse in 
the middle of the sequence to refocus low-frequency fluc- 
tuations. For Gaussian-distributed noise with a l//-type 
spectrum, the decay of both the FID and the Hahn-echo 
has the form p{t) = e~^/^^i ^-(t/T^^F/E)) ^ f^^st mea- 
suring Ti we can extract the pure dephasing times T^p 
andT^E [3]. 

In Fig. [sj^b) we plot T^^p and T^^e as a function of fz for 
the two values of A. We make a few observations: (i) The 
time T^E is 4-5 times longer than T^^f, consistent with the 
1 / /-type spectrum [25] . (ii) Both T^p and T^e fall off as 
|/z| is increased; this is because dEoi/dfz is zero only at 
/z = and increases approximately linearly with | fz \ over 
this range of fz. However, contrary to Refs. [3l[26], the 
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FIG. 3: (a) Energy decay time Ti. The decay does not 
depend strongly on the flux detuning, (b) Echo and free- 
induction decay times versus flux detuning. The position of 
the longest decay times is shifted to positive flux detuning, (c) 
Free-induction decay measured at /x = ±0.39. Depending on 
the sign of /x, the position of the longest decay time shifts in 
different directions. All measurements were done by sweeping 
fz while keeping /x constant. 



echo dephasing times are considerably shorter than the 
limit set by energy relaxation (2Ti) even close to /z = 0. 
This is a result of our device having an extra loop that 
controls A: fluctuations in /x will couple to A, which will 
couple to £^01 even when = 0. (iii) Close to /z = 0, 
the dephasing times are longer for A = 2.5 GHz than for 
A = 3.5 GHz: this is because the sensitivity dA/df^ to 
noise in /x is stronger at A = 3.5 GHz [see Fig.|2|c)]. (iv) 
The longest dephasing times do not occur at /z = but 
are slightly shifted to positive /z. The effect is visible for 
both values of A, but the shift is larger for A = 3.5 GHz. 

We investigate the shift in more detail by measuring 
the FID at A = 3.5 GHz for positive and negative values 
of /x. The results are shown in Fig.[3];c); at /x = -0.39, 
the longest dephasing times are shifted towards negative 
fz. Even though A is the same at both /x = 0.39 and 
/x = —0.39, the sensitivity dA/df^ is negative for /x < 
and positive for /x > 0. The fact that the shift depends 
on the sign of dA/df^ leads us suspect that there are 
correlations between /x and /z. We have checked that 
the shift does not depend on SQUID bias current for 
small excursions from the operating point [27j. 

The phase decay rate T^^ = is due to a 

combination of noise in both /x and /z. To calculate 
F^^F in the presence of correlations, we assume that the 
fluctuations Sf^ and 5fz are described by noise spec- 
tra of the form Sf^{uj) = (Ax^o^)/l^l ^/z(^) = 
{A^^q'^)/\(jj\. In addition we introduce the correla- 
tion spectrum Sf^f^{uj) = (Azx^o^)/l^l EHl EH] and 



(a) 20 

15 

CO 

10 

Ll_ 

5 




J-100 
^^-200 




1 V 1 V 1 










-aEo,/a/x- 


/^min(^zx 1) 


-5Eoi/afz- 



-0.5 

fz [xlO^] 



0.5 




FIG. 4: (a) Phase decay rate Fc^f expected from the sensitiv- 
ities shown in (b), with — — (3/i$o)- (b) Sensitivity 
of the qubit energy to flux fluctuations. For perfectly cor- 
related noise (czx = il)? the fluctuations cancel out when 
dEoi/df^ — ^dEoi/ dfz- (c) Measured decay rate Fc^f- The 
fit gives a correlation factor of Czx = —0.25. (d) Schematic 
picture of the surface-spin model. Spins located on the shared 
line marked by dc induce fields pointing in opposite direction 
in the two loops. 



the correlation coefficient Czx = ^/x) / \/ (^/z ) (^/x ) 
^zx/V^x^z- The total decay rate is \2b\ 
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The low- frequency cut-off a;iow/27r = 1 Hz is imposed by 
the measurement protocol. Figure [Ij^a) shows the result 
of Eq. Q , plotted for A^ = A^ and three different values 
of the correlation coefficient Czx- Without any correla- 
tions (czx = 0), the decay rate increases with |/z|, and the 
value at /z = is set by the noise in f^. However, if we 
assume perfectly correlated (czx = 1) or anti-correlated 
noise (czx = — 1)? the point of minimal decay shifts to 
negative and positive values of /z, respectively. In addi- 
tion, the decay rate goes to zero in those points. 

To understand how correlations can lead to shifts of 
the dephasing times, we consider how the fluctuations 
in /x and /z couple to the energy splitting Eqi. Fig- 
wce^h) shows the sensitivities dE^i/df^ and dEoi/dfz 
for A = 3.5 GHz, extracted from the data in Fig. [2j 
If the noise is completely correlated, the fluctuations in 
£^01 will perfectly cancel when dEoi/df^ = —dEoi/dfz. 
This happens at the point marked by Fmin(czx = 1) in 
Fig. ^h). On the other hand, if the noise is perfectly 
anti-correlated, the fluctuations in Eqi will cancel when 
dEoi/dfx = dEoi/dfz. This also explains why the phase 
decay data taken at /x = —0.39 shifts in the opposite 
direction [Fig.jsjc)]. Here, the derivative dA/df^ is neg- 
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ative, which gives dEoi/df^ a minus sign and makes the 
points rniin(czx = l) and rmin(czx = — 1) trade places. 

We can fit the measured decay rate to Eq. ([3| in or- 
der to extract the flux-noise amplitude and the correla- 
tion coefficient in our sample. The results are shown in 
Fig.gc), giving A, = (3.9/i$o)', = (3.1/i^o)' and 
c = —0.25 ± 0.05. We stress that the same fitting pa- 
rameters are used for both values of A. The reason why 
the data for A = 2.5 GHz shows a less dramatic shift 
in /z is because dEoi/df^ is smaller at this value of A. 
The values for dE^i/df^ and dE^i/df^ used in the fits 
were extracted from Fig. |2j giving de/df^ = 1.37 THz, 
dA/df^ = 110 GHz for A = 2.5 GHz and de/df^ = 
1.28 THz, SA/a/x = 142 GHz for A = 3.5 GHz. 

It is important to point out that the correlations dis- 
cussed so far have been between fluctuations J/z, Sf^ in 
the effective fluxes defined by Eq. ([T]). The fluctuations 
and the correlations in the geometric fluxes fi , /2 can be 
found by inverting Eq. ([T]). Assuming the fluctuations 
^/i, are described by spectra Sf-^^{uj) = (74i<I>q ^)/|cj| 
and Sf^{uj) = {A2%^)/\uj\, we find [39] 
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Ax = (3.1/i$o)', 



-0.55. 



(4) 



Fabrication imperfections may cause deviations in the 
sizes of the two junctions in loop 2, which will affect 
Eq. (16). From the results of Fig. [2|b), we know that 
the junction asymmetry is less than ±10%, which will 
lead to an uncertainty of the correlation coefficient Ci2 of 
about ±0.1. Note that the measurement is sensitive to 
correlations in the frequency range 1 Hz to 1 MHz, as 
set by the FID weighting function [25l HB] . 

The results of Eq. (16) show that the flux fluctuations 



in the two loops of the sample are anti-correlated. This 
is clear evidence that the noise is not due to a global fluc- 
tuating magnetic field, since such fluctuations would give 
positive correlations. Instead, we consider a model where 
the noise is due to a large number of randomly oriented 
spins distributed over the metal surface [see Fig. |4]^d)]. 
The spins will generate local magnetic fields that is picked 
up by the loops, and fluctuations of the spin orientations 
will give rise to flux noise [7 . In our sample, spins that 
are located on the line that is shared by the loops [marked 
by dc in Fig.jlj^d)] will induce fields of opposite directions 



in loop 1 and loop 2. This will give rise to negative cor- 
relations between the fluxes in the loops. 

To get an estimate of the expected correlations, we 
assume (i) that the field of a single spin only couples 
to a loop if it is sitting on that loop's circumference, 

(ii) that the width of the sample wires is constant, and 

(iii) that the amount of flux that a loop picks up from 
a single spin is independent of the spin's position along 
the width of the wire. The assumptions seem reasonable 
based on numerical simulations of the coupling between 
a SQUID and the magnetic moment of a single spin [7]. 
Under these conditions the flux fluctuations generated 
by an ensemble of randomly oriented spins will scale as 
5f oc yfpd^ where p is the spin density and d the length 
of the loop. The correlation coefficient becomes 
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Here dc is the length of the line segment shared by both 
loops, and di^ ^2 are the lengths of remaining sections of 
the two loops [Fig. [4]^d)]. With values relevant for our 
sample {di = 10.8 /im, d2 = 6.0 /im, dc = 4.8 /im), we 
find Ci2 = —0.37. This is smaller than the value extracted 
experimentally, but the agreement is reasonable consid- 
ering the simplicity of the model. Further work is needed 
to develop a more realistic model of the device [24 , and 
to describe the correlations when spin interactions are 
taken into account [151 US] • 

To conclude, we have investigated flux noise corre- 
lations by measuring dephasing rates in a tunable flux 
qubit. We find the flux fluctuations in neighboring loops 
to be anti-correlated, in agreement with predictions from 
models where the flux noise is generated by randomly ori- 
ented spins on the metal surface. The ability to extract 
noise correlations provides important information about 
the microscopic origin of the flux noise. The method pre- 
sented here is general and can be extended to samples 
with different geometries. Combining this method with 
recent improvements in multi-pulse noise spectroscopy 
[26 opens up the possibility to measure correlations at 
different frequencies and thereby probe the dynamics of 
the surface spin system. 

We thank T. Orlando and R. McDermott for helpful 
discussions, K. Harrabi for assistance with device fabri- 
cation and P. Forn-Diaz for performing numerical simu- 
lations of the qubit spectrum. 
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Supplementary material 

Device drawing 

The device was fabricated using standard two-angle 
shadow evaporation techniques. The two layers of metal 
are deposited from different angles, resulting in overlap- 
ping superconductors separated by a thin film of insu- 
lating aluminum oxide. The method is used to form the 
Josephson junctions in the system (see Fig. [5|. 




FIG. 5: Drawing of the device, showing the location of metal 
from the two evaporation stages. 



SEM picture of the device 

Figure [6] shows a SEM picture of a device with the same 
geometry as the one used in the experiments. The qubit 
and the read-out SQUID are in the upper-right part of 
the figure. The two lines next to the qubit are the local 
flux bias lines, while the thicker line in the lower-left part 
of the figure is the microwave drive line. 



Derivation of flux correlations 

We start by defining the normalized fluxes /i = ^i/^o 
and /2 = $2/^0 in the two qubit loops. Correlations 
between fluctuations Sfi and 6/2 in those parameters are 
characterized by the correlation coefficient C12, defined as 



C12 



The qubit parameters e = e{fz) and A 
controlled by the effective fluxes 

/z = (/i + /2/2) - 1/2, 

/x = /2- 



A(/x) 



(6) 



are 
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FIG. 6: SEM picture of a device with the same geometry as 
the one used in the experiments. 



The coefficient for describing correlations between the 
fluctuations Sf^ and 6f^ in the effective fluxes is given by 



{SflSf2) 



VW)W) 



(8) 



We want to be able to extract C12 from the measured 
correlations Czx between the effective fluxes. We start 
by inserting the expression for {8fi 6/2) from Eq. (|6| into 
Eq. (Is]). This gives 




(9) 



Solving for C12 gives 



C12 = Cz 



'Ml 



m 

(10) 

In addition, we need to express the fluctuations {dfl) in 
terms of the effective fluctuations Czx, (^/x) (^/z)- 
From the definition in Eq. ([7|), we have 

(O = mi+sf2/2f) 

Using the expression for 



W!)+{ShSf2) + {6fl)/4. (11) 
fi ^12) from Eq. (|6|, we get 

^m) m) + (<^/i>/4. (12) 



Inserting the expression for C12 from Eq. ( |1Q[ ) gives 



(13) 
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Noting that (<5/|) = (<5/2), we have 



w!) = (sf!) + c^^vrnm- (c')/4. (14) 

Finally, solving for {5fi) gives 



(<5/i ) = (O + (<5/x')/4 - c^.vmm- (15) 

We assume that the fluctuations J/i, (5/2, % and J/x 
are all described by noise spectra of the form Sf^{uo) = 
A,/\co\, Sf.iio) = A2/\co\, SM = and SM = 



Ax/|w|. Using this and summarizing the results of 
Eq. ([lO| and Eq. we get 



A2 = ^x, 

Az 1 / 
aI~ 2\Ia'i' 



C12 



This is Eq. (4) in the main paper. 
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